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ABSTRACT 

We present the first phase-coherent measurement of a braking index for the 
young, energetic rotation-powered pulsar PSR J1846— 0258. This 324ms pulsar 
is located at the center of the supernova remnant Kes 75 and has a characteristic 
age of t c = 723 years, a spin-down energy of E = 8.3 x 10 36 ergs _1 , and inferred 
magnetic field of 4.9 x 10 13 G. Two independent phase-coherent timing solutions 
are derived which together span 5.5 yr of data obtained with the Rossi X-ray 
Timing Explorer. In addition, a partially phase-coherent timing analysis con- 
firms the fully phase-coherent result. The measured value of the braking index, 
n = 2.65 ± 0.01, is significantly less than 3, the value expected from magnetic 
dipole radiation, implying another physical process must contribute to the pul- 
sar's rotational evolution. Assuming the braking index has been constant since 
birth, we place an upper limit on the spin-down age of PSR J1846— 0258 of 884 yr, 
the smallest age estimate of any rotation-powered pulsar. 

Subject headings: pulsars: general — pulsars: individual (PSR J1846— 0258) — X- 
rays: stars 
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1. Introduction 

The measurement of pulsar braking indices (n) is crucial to the understanding of the 
physics underlying pulsar spin down, assumed to be of the form 

v = -Kv n , (1) 

where v is the pulse frequency, v is the frequency derivative and K, assumed to be constant, 
is related to the pulsar's magnetic field and moment of inertia. By taking a time derivative 
of the above spin-down equation, we find that n = ui)/u 2 , where i) is the second frequency 
derivative. It is typically assumed that pulsars radiate as perfect magnetic dipoles, implying 
n = 3. This assumption is used implicitly for the estimation of pulsar magnetic fields as well 
as to calculate characteristic ages (defined as r c = P/2P corresponding to n — 3). However, 
of the five unambiguous measurements of pulsar braking indices obtained so far, all yield 
a value of n < 3 (Lyne et al. 1993, 1996; Camilo et al. 2000; Livingstone et al. 2005a,b). 
Explanations for this discrepancy include the possibility that: the relativistic pulsar wind 
affects the spin-down (Michel & Tucker 1969); the pulsar may suffer a propeller torque from a 
putative supernova fallback disk (Menou et al. 2001); the pulsar has a time- varying magnetic 
moment (Blandford & Romani 1988); or, the pulsar cannot be modelled as a point dipole, 
rather, a dipole of some finite size which leads to n < 3 (Melatos 1997). 

There are few pulsars that are potential candidates for a significant measurement of n; 
4 of the 5 measured values are from sources with characteristic ages under 2000 yr. The 
first requirement for a significant measurement of n is that the pulsar spins down sufficiently 
quickly for a measurement of i). The second requirement is that the position of the pulsar 
is accurately known at the ~1 arcsecond level. The third requirement is that the spin- 
down must not be seriously affected by glitches, sudden spin-ups of the pulsar, or timing 
noise, a long-term, low-frequency stochastic wandering of the rotation about the overall 
trend. Typically, glitches begin to seriously affect smooth spin down at characteristic ages 
of ~ 5 — lOkyr (McKenna & Lyne 1990; Marshall et al. 2004). Thus many of the pulsars 
that may spin down fast enough for a measurement of n are irretrievably contaminated by 
glitches (e.g. Marshall et al. 2004). Timing noise varies from object to object, though is 
roughly correlated with spin-down rate (e.g. Arzoumanian et al. 1994) and can prevent a 
measurement of n in a finite data set in an unpredictable way. 

The very young, energetic pulsar PSR J1846— 0258 was discovered at the center of the 
supernova remnant Kes 75 (Gotthelf et al. 2000). Neutral hydrogen absorption measurements 
indicate that the pulsar resides well across the Galaxy, roughly at a distance of 19 kpc (Becker 
& Helfand 1984). Radiation from this pulsar is detected exclusively in the X-ray band, with a 
0.5-10 keV luminosity of L x = 4.1 x 10 35 ergs -1 for a distance of 19 kpc (Helfand et al. 2003). 
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If correct, the spin-down energy to X-ray conversion efficiency is 1.6%, 6 times greater than 
that found in the Crab pulsar. The X-ray spectrum is typical of rotation-powered pulsars, 
with non-thermal power-law emission with T = 1.4 (Helfand et al. 2003). However, its large 
inferred magnetic field of B = 4.9 x 10 13 G and spectrum place it in an emerging class 
of rotation-powered pulsars with magnetar-strength fields (e.g. Kaspi & McLaughlin 2005). 
With a characteristic age of r c = 723 yr, PSR J1846— 0258 is likely the youngest of all known 
rotation-powered pulsars. 

In this paper, we present the first phase-coherent timing solution for the Kes 75 pulsar 
based on 5.5 yr of X-ray timing data from a long-term monitoring campaign with the Rossi 
X-ray Timing Explorer (RXTE). 

2. RXTE Observations and Analysis 

Observations of PSR J1846— 0258 were made using the Proportional Counter Array 
(PCA; Jahoda et al. 1996) on board RXTE. The PCA consists of an array of five col- 
limated xenon/methane multi-anode proportional counter units (PCUs) operating in the 
2-60 keV range, with a total effective area of approximately 6500 cm 2 and a field of view 
of ~ 1° FWHM. 

The data span 6.3 yr from 1999 April 18 - 2005 July 27 and were collected in "GoodXenon" 
mode, which records the arrival time (with l-/is resolution) and energy (256 channel resolu- 
tion) of every unrejected event. Typically, 3 PCUs were operational during any observation. 
We used all layers of each operational PCU in the 2-60 keV range, as this maximizes the 
signal-to-noise ratio for this source. We analyzed 171 observations, resulting in 158 detec- 
tions of the pulse for a total integration time of 337 hr. The data were unevenly spaced 
throughout the 6.3 yr of observations, as shown in Figure 1. Because PSR J1846— 0258 was 
occasionally not the primary target of RXTE in the observations we have used, integration 
times ranged from ~1.5 to >25 ks, resulting in a variety of signal-to-noise ratios for individual 
pulse profiles. 

Observations were downloaded from the HEASARC archive 1 and data from the different 
PCUs were merged and binned at (1/1024) ms resolution using the Ftools 2 'seextrct' and 
'fselect'. Photon arrival times were converted to barycentric dynamical time (TDB) at 
the solar system barycenter using the J2000 source position RA = 18 h 46 m 24!94, Dec = 



^ttpi/yheasarc. gsfc.nasa.gov/docs/archive. html 
2 http://heasarc. gsfc.nasa.gov/docs/software/ftools/ 
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-02°58'30.1" (Helfand et al. 2003) and the JPL DE200 solar system ephemeris with the FITS 
tool 'faxbary'. The remainder of the analysis was performed using software developed at MIT 
for handling RXTE data. We determined an initial ephemeris by merging adjacent short 
observations and performing a periodogram analysis that resulted in a set of 134 frequency 
measurements spanning the 6.3-yr monitoring interval. From these measurements, we derived 
an initial frequency derivative, v. This ephemeris was used to fold each time series with 16 
phase bins. This number of bins was chosen because of the lack of features in the roughly 
sinusoidal profile and the resulting reasonable signal-to-noise ratio for individual profiles, 
particularly those made from short integrations. Resulting profiles were cross-correlated in 
the Fourier domain with a high signal-to-noise ratio template created by adding phase- aligned 
profiles from all observations, shown in Figure 2. The cross- correlation process assumes that 
the pulse profile is stable; indeed, we found no evidence for variability. We implemented 
a Fourier domain filter by using only the first six harmonics in the cross-correlation. For 
each observation, the cross-correlation yeilded the time of arrival (TOA) of phase- zero of the 
average pulse profile at the fold epoch. The TOAs were fitted to a timing model (see §2.1 
and §2.2) using the pulsar timing software package TEMPO 3 . After phase-connecting the 
data (see §2.1), we merged observations occurring on a single day and used the ephemeris 
to re-fold the data in order to obtain more accurate TOAs. This process produced 81 TOAs 
with a typical uncertainty ~8ms (~2.5% of the pulse period), improved from a typical TOA 
uncertainty of ~10 ms. 



2.1. Phase-Coherent Timing Analysis 

To obtain a phase- coherent timing solution, each turn of the pulsar is accounted for by 
fitting TOAs to a Taylor expansion of pulse phase, 0. At time t, <fi can be expressed as: 

<P(t) = 0(t o ) + v (t - t ) + \vv{t - t ) 2 + -i> (t - t ) 3 + ..., (2) 

where the subscript denotes a parameter evaluated at the reference epoch, t . TOAs and 
initial parameters are input into TEMPO, which gives as output refined spin parameters 
and residuals. 

To determine the spin parameters for PSR J1846— 0258, we obtained phase-connected 
timing solutions for 5.5 yr of X-ray timing data spanning 2000 January 31 to 2005 July 27, 
and including 78 TOAs. The sampling of the observations over 6.3 yr includes several large 
gaps and that precludes a single phase-coherent timing solution for the entire interval. 



3 http: / / www.atnf.csiro.au/research/pulsar / tempo 
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The first RXTE observations of the pulsar (which resulted in its discovery) occurred 
~9 months prior to the commencement of regular monitoring observations, and so were not 
useful in our phase-coherent analysis. The first coherent solution is valid over MJD 51574- 
52837 (3.5 yr), while the second coherent solution is valid over MJD 52915-53578 (1.8 yr), as 
indicated in Figure 1. 

We used our initial ephemeris (described in the previous section) to bootstrap a phase- 
coherent solution, valid over the 3.5 yr interval from MJD 51574-52837. This solution includes 
z/, z> and i), whose values are given in Table 1. In the process of phase connection, it became 
clear that a small glitch occurred at MJD 52210±10. Our measured glitch parameters are 
Ais/v = 2.5(2) x 1CT 9 and Az>/z> ~ 9.3(1) x 10~ 4 , as determined with the glitch fitting facility 
in TEMPO. The relatively wide spacing of data near the glitch epoch prevent the detection 
of any glitch recovery. In fact, it is possible that the initial frequency jump was larger and 
recovered significantly before the following observation. Timing residuals after subtraction 
of our best-fit timing model, including the glitch, are shown in the top panel of Figure 3. 
Note that systematic trends remaining in the residuals are likely the result of timing noise, 
common among young pulsars, though unmodelled glitch recovery may also contribute to 
the observed residuals. Timing noise processes are known to contaminate measured spin 
parameters, hence it is typically advisable to remove the systematics from the residuals by 
fitting additional frequency derivatives until the residuals are consistent with Gaussian noise 
(e.g. Kaspi et al. 1994). For this pulsar, a total of eight frequency derivatives were required 
to obtain Gaussian distributed residuals. Timing residuals with all eight derivatives removed 
are shown in the bottom panel of Figure 3. Fitting additional frequency derivatives improves 
the x 2 from 2933 for 43 degrees of freedom to 77 for 37 degrees of freedom. This x 2 value 
indicates that the fit does not completely describe the data, however, this is not uncommon 
when fitting timing noise, which is often not well described by a polynomial. The braking 
index, resulting from this 'whitened' timing solution, is n = 2.64 ±0.01. Deterministic spin- 
down parameters (i.e. not the higher order derivatives that represent timing noise) as well 
as glitch parameters for this timing solution are given in Table 1. 

Phase was lost over a 78-day gap in the data beginning at MJD 52837, made clear by the 
fact that a solution attempting to connect over this gap fails to predict the pulse frequency 
at previous epochs. This loss of phase could be due either to timing noise or another glitch. 
However, the estimated change in frequency over the gap, calculated from phase-coherent 
solutions on each side of the gap, is negative, ie. in the opposite direction from a conventional 
glitch. This implies that a glitch alone cannot account for the loss in phase. 

A second phase-coherent solution was obtained for the 1.8-yr interval from MJD 52915- 
53579, with z/, z>, i) fitted. The residuals after subtraction of the best-fit parameters are 
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shown in the top panel of Figure 4. As systematic trends, again interpreted as timing 
noise or possibly unmodelled glitch recovery, remained in the residuals, three higher-order 
frequency derivatives were fitted. The resulting residuals are shown in the bottom panel 
of Figure 4. The measured braking index resulting from the 'whitened' timing solution is 
n = 2.68 ± 0.03, in agreement with that measured from the first segment of timing data. 
Deterministic parameters for the second timing solution are given in Table 1. 

2.2. Partially Phase-Coherent Timing Analysis 

In order to mitigate the effects of timing noise, we also performed a partially phase 
coherent analysis. In this way we obtained measurements of v spanning 6.3 yr of data as 
well as v and i) over 5.5 yr of data. This method is useful to detect small glitches, as well 
as to obtain more accurate measurements of n in some cases (e.g. Livingstone et al. 2005a). 
Using the overall phase- coherent ephemeris as a starting point, closely spaced observations 
were phase-connected to obtain a local measurement of v. In this way we obtained a total of 
22 measurements. A two degree polynomial was fitted to these measurements to get another 
measurement of n, however, in this case, the glitch near MJD 52210 and the possible glitch 
between MJD 52837-52915 seriously restricted the available time baseline and rendered this 
analysis of limited value. The most constraining measurement is n = 2.83 ± 0.39, from 13 
values of v spanning MJD 51286-52199, in agreement with our phase-coherent value. 

We repeated this process, measuring eight independent values of z>, shown in the top 
panel of Figure 5. Note that the first unambiguous z> measurement occurred before the 
glitch near MJD 52210 and the glitch can clearly be seen in the Figure. The bottom panel 
of Figure 5 shows the post-glitch slope removed from the data, highlighting the change in v 
at the time of the glitch, Ai>/z> = (9.5 ± 0.3) x 10~ 4 , in agreement with the value obtained 
from the phase-coherent fit. 

Note that at the 78 day data gap between MJD 52837 and 52915, there is no clear Az> 
visible in Figure 5. This indicates that if a glitch did occur during this period, it consisted 
of primarily a change in v. The only glitch consistent with the data would have had to have 
occurred near MJD 52910 with magnitude Az//V < 5 x 10~ 8 . A glitch with no change in Az> 
is typical of pulsars with larger characteristic ages (> 5kyr), but cannot be ruled out in this 
case. 

We also performed two linear weighted least square fits on the z> measurements, both 
before and after the gap at MJD 52837-52915, to obtain measurements of v and thus n. The 
first segment resulted in n = 2.47 ± 0.07, while the second resulted in n = 2.65 ± 0.14, 1.3a 
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from each other. Uncertainties were obtained from a bootstrap analysis, since we suspected 
that the formal uncertainty may underestimate the true uncertainty because of the presence 
of timing noise. The bootstrap is a robust method of determining errors when a small number 
of sample points is available, as in this case (Efron 1979). The first measurement is 2.5a 
from our phase-coherent average value of n — 2.65 ± 0.01, while the second measurement is 
in agreement with it. 

Finally, we obtained four independent measurements of i) by phase-connecting as much 
data as possible with just u, v and is, such that there was no noticeable contamination by 
timing noise appearing in the resulting residuals. The resulting n measurements are shown 
in Figure 6. Note that there is some indication that n is growing with time, at at rate 
of An/ At ~ 0.02 ± 0.01 yr -1 , resulting in a change in n of ~ 5% over 5.5 yr. However, 
this effect is only significant at the 2a level. A similar analysis of 21 yr of data of the 
young pulsar PSR B1509— 58 showed long-term variations of ~1.5% in n, though these were 
shown to fluctuate about an average value of n (Livingstone et al. 2005b). Further timing 
observations of PSR J1846— 0258 will show if its value of n is changing secularly or is subject 
to timing noise variations as in PSR B1509— 58. 

3. ASCA and BeppoSAX Observations 

In order to check the validity of our RXTE analysis, we re-analyzed the timing signal in 
two archival observations of the pulsar obtained with the Advanced Satellite for Cosmology 
and Astrophysics (ASCA) observatory (Tanaka et al. 1994). Details of the first observation 
(1993 October) can be found in Blanton & Helfand (1996) and Gotthelf et al. (2000), while 
details of the second observation (1999 March) are described in Vasisht et al. (2000). Using 
the prescription in the ASCA Data Reduction Guide 4 and standard FITS tools, we extracted 
photons from the two gas-imaging spectrometers (GISs) in the 3-8 keV range from a 4' radius 
aperture surrounding the source. We chose this energy range in order to preferentially select 
the harder energy photons which mainly come from the pulsar rather than the softer energy 
events having a high probability coming from the supernova remnant (Gotthelf et al. 2000). 
We adjusted photon arrival times to the solar system barycenter and then performed a 
periodogram analysis on the event data centered on the pulse frequency predicted by our 
RXTE ephemeris. We folded the data with 10 phase bins and detected a weak pulse in each 
observation. We generated a x 2 periodogram in frequency space and identified the peak 
v. To determine the uncertainty in the measured frequencies, we performed a Monte Carlo 



4 http://heasarc. gsfc.nasa.gov/docs/asca/abc/abc. html 
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simulation that created fake data sets containing a periodic signal (at the v and amplitude 
determined from the periodogram) and Poisson noise. A periodogram was then performed on 
each noisy fake data set and the best v was recorded. We performed 500 iterations for each 
observation and took the standard deviation as the la uncertainty on our v measurement. 
Using this method, we determined that for the 1993 October observation, v = 3.094516(2) Hz 
and for the 1999 March observation, v = 3.082859(3) Hz. 

To verify our RXTE results, we compared the above ASCA- measured frequencies, as 
well as two 1999 BeppoSAX timing observations reported on by Mereghetti et al. (2002), with 
the RXTE ephemeris prediction. Both BeppoSAX v measurements agree with our ephemeris 
within uncertainties, as do the two ASCA frequencies. That the 1993 ASCA measurement 
of v agrees with our prediction is somewhat surprising, since the ephemeris used to predict 
the value covers dates beginning seven years later. This does not necessarily exclude the 
possibility of glitches having occured in the interim however, since glitches in other young 
pulsars have been measured that sometimes have amplitudes ~ 3 orders of magnitude smaller 
than we could detect given the uncertainties here (Livingstone et al. 2005a). 

If no glitches occurred during the 6-yr gap from 1993 to 1999, then it is possible to fit 
all frequency measurements prior to the MJD 52210 glitch with a second order polynomial 
in order to measure the braking index. Performing this analysis, the resulting braking index 
is n = 2.44 ± 0.09, 2.2<r from our phase-coherent result. 



4. Discussion 

Establishing true ages of pulsars is important for several aspects of neutron-star studies, 
including neutron star cooling, population synthesis, spatial velocity estimates, and to con- 
sider their associations with supernova remnants. Based on spin-down properties alone, the 
true pulsar age is impossible to determine. However, we can estimate the age based on the 
standard spin-down model (z> oc — v n ) given our measurement of n. Integrating this model 
yields, 

1 v I v 



Assuming that the initial spin frequency was much larger than the current value, u Q ^> u, 
the age estimate for PSR J1846— 0258 becomes, 

r< — -< 884yr. (4) 

~ n- lz> ~ J w 

Since the initial spin frequency is unknown, the above estimate provides an upper limit on 

the pulsar age. Therefore, assuming that n and K have remained constant over the lifetime 
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of the pulsar, the upper limit on the age of PSR J1846— 0258 is ~884yr. This value is larger 
than the characteristic age of 723 yr, but still smaller than the unambiguously known age of 
the next youngest pulsar, that in the Crab Nebula. 

Mereghetti et al. (2002) did an analysis of frequency measurements of the ASCA and 
BeppoSAX data, as well as a small subset of the RXTE data reported here. Although 
data gaps prevented them from a firm conclusion regarding n because of the possibility of 
glitches, their tentative estimate (made assuming no glitches nor substantial timing noise) 
was n = 1.86 ± 0.08, inconsistent with our result. There are two plausible explanations for 
this discrepancy. The first is that significant timing noise and/or glitches cannot be definitely 
excluded in the 6-yr gap between the successive ASCA observations. The second is that the 
coordinate used to barycenter these data was ultimately found to be off by ~ 7". The earlier 
study used coordinates derived from ASCA data with a 20" position uncertainty, prior to 
the availability of the new measurement provided by the Chandra X-ray Observatory used 
herein. 

Our measured n = 2.65±0.01 is significantly less than 3, as is the case for all established 
values of n, which are shown in Table 2. There are several ideas for the nature of the deviation 
from the prediction of simple magnetic dipole braking. A time-varying magnetic moment 
could produce an observed n less than the true Uq = 3 (Blandford 1994). A varying i?-field 
can in principle be verified with a measurement of the third frequency derivative, is, which 
may or may not ultimately be possible for this source, depending on the number and strength 
of the glitches it experiences, as well as on the strength of its timing noise. Presently, v 
has been measured for only two pulsars, both of which are consistent with a constant B 
(Lyne et al. 1993; Livingstone et al. 2005b). Another suggested explanation for n < 3 is 
that a fallback disk forms from supernova material and modulates the spin-down of young 
pulsars via a propeller torque. Spin-down via a combination of magnetic dipole radiation 
and propeller torque results in 2 < n < 3 (e.g. Alpar et al. 2001). However, this requires that 
the disk material does not suppress the pulsed emission during the propeller phase, which 
is difficult to achieve (Menou et al. 2001). Angular momentum loss due to a stellar wind 
would result in n — 1 (Michel & Tucker 1969). Thus spin-down due to some combination 
of relativistic pulsar winds associated with young neutron stars and observed indirectly as 
pulsar wind nebulae (e.g. Roberts et al. 2003), and magnetic dipole spin down may result in 
measured values n < 3. 

Melatos (1997) presented an intriguing solution to the n < 3 problem. He postulated 
that the radius pertinent to dipole radiation is not the physical neutron-star radius, but a 
"vacuum radius" , associated with the location closest to the neutron star where field-aligned 
flow breaks down. Since this radius can be significantly larger than the neutron-star radius, 
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the system can no longer be treated as a point dipole and n = 3 is not necessarily true. His 
model provides a prediction for n, given three observables: u, v and a, the angle between 
the spin and magnetic axis. This model predicts that 2 < n < 3 for all pulsars and that n 
approaches 3 as a pulsar ages. Currently, there is no measurement of a for PSR J1846— 0258. 
However, for the model to be consistent with our measured n within 3a, a must lie between 
8.1° and 9.6°. This is small but not unreasonable given the very broad pulse profile for 
this source. Future radio polarimetric observations could in principle constrain a, though at 
present, no radio detection of this source has been reported (Kaspi et al. 1996). 

The measurement of n for PSR J1846— 0258 brings the total number of measured brak- 
ing indices to six, shown in Table 2, along with v, z>, r, r c , -Bdipoic and E for each object. 
Comparing this new value to the other four measurements obtained via typical timing meth- 
ods (i.e. excluding the Vela pulsar whose n was measured using a different method, due to 
the large glitches experienced by this pulsar), we find no correlations between n and any of 
the other parameters. If the Vela pulsar is included in the analysis, there is a slight corre- 
lation between n and characteristic age. However, there is a relatively large scatter among 
the five younger pulsars, suggesting the Vela pulsar's value could also be just an extremum 
of the scatter. 

The large scatter in the observed values of n could be intrinsic to the pulsars, for 
instance, there may exist a relationship between n and a or B, as suggested by Melatos 
(1997). One way to test this possibility would be to obtain more precise measurements of 
a for the youngest pulsars, e.g. via polarimetric observations. This, however, may prove 
difficult given that young radio pulsars tend to show very flat position angle swings (e.g. 
Crawford et al. 2001; Johnston et al. 2005). Alternatively, the scatter could be the result 
of a physical process outside the neutron star that is affecting pulsar spin down, such as a 
supernova fallback disk as suggested by Menou et al. (2001). A supernova fallback disk was 
recently discovered by Wang et al. (2006) around a young neutron star, though there is no 
indication that the disk and the neutron star are interacting. 

PSR J1846— 0258 is in an emerging class of high magnetic field rotation-powered pulsars 
(e.g. Camilo et al. 2000; Pivovaroff et al. 2000; Kaspi & McLaughlin 2005; McLaughlin et al. 
2003). The existence of these sources raises the question of why they appear to be rotation 
powered instead of powered by the decay of their large magnetic fields, as are the magnetars. 
The magnetar model proposed by Thompson et al. (2002) provides a prediction for the X-ray 
luminosity of a neutron star with a magnetic field of ~ 10 14 G given a measurement of n. In 
the model, the neutron-star magnetosphere suffers a large scale 'twist' of the north magnetic 
hemisphere with respect to the southern hemisphere, and resulting magnetospheric currents 
both scatter thermal surface photons and themselves heat the surface upon impact. Both 
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effects result in the observed X-ray emission. The model predicts a one-to-one relationship 
between twist angle and n, as well as L x . In the case of PSR J1846— 0258, if the above 
model were applicable, the measurement of n — 2.65 + / — 0.01 would imply a twist angle of 
<f> ~ 1.2 rad, which would predict an X-ray luminosity of L x ~ 3.7x 10 35 ergs -1 . The observed 
luminosity of PSR J1846— 0258 in the 0.5-10keV energy range is L x = 4.1 x 10 35 ergss _1 
for a distance of 19kpc (Helfand et al. 2003). This value is consistent with the magnetar 
model prediction. This suggests that the pulsar may actually be a magnetar. However, 
the estimated distance is is uncertain by a factor of two and likely overestimated, since the 
remnant flux and size are much larger than expected for the pulsar's age. Thus the true 
X-ray luminosity of PSR J1846— 0258 could be significantly less than the magnetar model 
prediction, though this will likely remain uncertain. In any event, since the observed X-ray 
flux can be easily accounted for by the spin-down luminosity of the pulsar, PSR J1846— 0258 
is clearly not an 'anomalous' X-ray pulsar. Moreover, the pulsar's X-ray spectrum is much 
harder than that of any anomalous X-ray pulsar, and inconsistent with what is predicted in 
the magnetar model. Indeed the spectrum is much more in line with those seen from rotation- 
powered pulsars (Gotthelf 2003). More theoretical work needs to be done to explain the 
difference between the magnetars and the rotation-powered pulsars having inferred magnetar- 
strength fields. 
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Fig. 1.— Distribution of RXTE observations of PSR J1846-0258 over the 6.3 yr interval. 
Each bar represents all observations occurring on a single day. 81 merged observations are 
shown in total. The mid-length bars represent data included in the first phase-connected 
solution, while the short bars represent data included in the second phase-connected solution. 
The initial RXTE observations (tallest bars) of the source occur significantly before the 
monitoring observations and cannot be unambiguously phase connected with the later data. 
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Fig. 2. — Integrated pulse profile of 337 hr of data over 6.3-yr of RXTE observations of 
PSR J1846— 0258. Two full pulse cycles are shown for clarity. 
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Fig. 3. — Phase-coherent X-ray timing analysis of the young pulsar PSR J1846— 0258 span- 
ning a 3.5-yr interval from MJD 51574 -52837. Top panel: Residuals with u, i>, i), as well 
as glitch parameters Au and Az> fitted. The glitch epoch, MJD 52210 is indicated by the 
arrow. Bottom panel: Residuals with glitch parameters and eight frequency derivatives in 
total fitted to render the residuals consistent with Gaussian noise. Fitting the additional 
parameters improves the x 2 from 2933 for 43 degrees of freedom to 77 for 37 degrees of 
freedom. Although these parameters do not completely describe the data, the x 2 is not 
improved significantly by fitting additional frequency derivatives, which is not uncommon 
when timing noise is present. 
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Fig. 4. — Phase-coherent X-ray timing analysis of PSR J1846— 0258 spanning a 1.8-yr interval 
from MJD 52915-53579. Top panel: Residuals with u, z>, and v fitted. Bottom panel: 
Residuals with five frequency derivatives total fitted to render the residuals consistent with 
Gaussian noise. Fitting the additional parameters improves the \ 2 from 271 for 27 degrees 
of freedom to 26 for 24 degrees of freedom. 



-18- 



Epoch (Years) 
2002 2003 2004 



\ 

N 

O 

G 0) 
CD > 

^ ■+-> 

ST « 

CD k, 

PL, ^ 

CD 
Q 



-6.68X10" 11 - 



-6.69X10" 11 



CD 
> 

CO 



CD ' 

Q W 

g 2 

^ CO 

&< u 

CD 5-1 
S-i 

[in 



-6.7x10 



-li _ 



-6.71X10" 11 



8xl0" 14 - 



6xl0" 14 
4xl0" 14 b 
2xl0- 14 


-2xl0- 14 



i — i 1 1 1 1 — i 1 1 — i 1 1 — | 1 — ' , ^ ',~ 



i — # — i 

i — * — i 
H 1 1— 



I 3E 1 



J I L_ 



J_ 



I * 1 



I 1» 1 



I Z 1 I * I 



_l I I L_ 



i — s — i i— i— i 



52500 53000 
Epoch (MJD) 



_l I I L 



Fig. 5. — Eight phase-coherent i> measurements spanning 5.5 yr of RXTE monitoring obser- 
vations of PSR J1846— 0258. Top panel: Measurements of v. There is a clear discontinuity 
after the first measurement of z>, which we interpret as a glitch. Bottom panel: v measure- 
ments with the post-glitch slope fitted to highlight the discontinuity, which is quantified as 
Az>/z>= (9.5 ±0.3) x 10" 4 . 
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Fig. 6. — Phase-coherent measurements of the braking index, n of PSR J1846— 0258. Each 
phase-connected solution has only v, v and v fitted and has Gaussian residuals. The braking 
index is consistent with being constant, since the slight increasing trend of An/ At ~ 0.02 ± 
0.01 yr _1 is only significant of the 2a level. 
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Table I. Phase-coherent Timing Parameters for PSR J1846— 0258. 
First phase-coherent solution a 



Dates (Modified Julian Day) 

Dates (Years) 

Number of TOAs 

Epoch (Modified Julian Day) 

v (Hz) 

v (10- 11 s- 2 ) 

v (10~ 21 s" 3 ) 

Braking Index, n 

Number of derivatives fitted 

RMS residuals (ms) 

Glitch epoch (Modified Julian Day) 

Au/u 

Az>/z> 



51574.2 - 52837.4 
2000 Jan 31 - 2003 Jul 17 
47 
52064.0 
3.0782148166(9) 
-6.71563(1) 
3.87(2) 
2.64(1) 
8 
7.65 
52210(10) 
2.5(2) x 10~ 9 
9.3(1) x 10" 4 



Second phase-coherent solution 



Dates (Modified Julian Day) 

Dates (Years) 

Number of TOAs 

Epoch (Modified Julian Day) 

v (Hz) 

v (10- 11 s- 2 ) 

v (10- 21 s" 3 ) 

Braking Index, n 

Number of derivatives fitted 

RMS residuals (ms) 



52915.8 - 53578.6 
2003 Oct 3 - 2005 Jul 27 
31 
53404.0 
3.070458592(1) 
-6.67793(5) 
3.89(4) 
2.68(3) 
5 
5.61 



a Quoted uncertainties are formal la uncertainties as reported by 
TEMPO. 
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Table 2. Spin and inferred parameters for pulsars with measured n ordered by spin-down 
age 



Name 


n a 


V 


V 






R d 

-Ddipolc 


E e 


Ref 








(10- n s- 2 ) 


(yr) 


(yr) 


(10 12 G) 


(10 36 erg s- 1 ) 




J1846-0258 


2.65(1) 


3.07 


-6.68 


723 


884 


49 


8.1 


(1) 


B0531+21 


2.51(1) 


30.2 


-38.6 


1240 


1640 


3.8 


460 


(2) 


B1509-58 


2.839(3) 


6.63 


-6.76 


1550 


1690 


15 


18 


(3) 


J1119-6127 


2.91(5) 


2.45 


-2.42 


1610 


1680 


42 


2.3 


(4) 


B0540-69 


2.140(9) 


19.8 


-18.8 


1670 


2940 


5.1 


150 


(5) 


B0833-45 s 


1.4(2) 


11.2 


-1.57 


11300 


57000 


3.4 


6.9 


(6) 



a Uncertainties on n are in the last digit. 
Characteristic age is given by r c = v/2i>. 

inferred upper limit timing age given n, r = v / (n — l)z>, assuming z/j >> v. 

d Dipole magnetic field estimated by -Bdipoic = 3.2 x 10 19 (— u/u 3 ) G, assuming a = 90° and 
n — 3. 

c Spin-down luminosity, E = 47r 2 Jz/z>, where it is assumed that I = 10 45 gcm 2 for all pulsars. 

'References (1) This work, (2) Lyne et al. (1993), (3) Livingstone et al. (2005b), (4) Camilo 
et al. (2000), (5) Livingstone et al. (2005a) (6) Lyne et al. (1996). 

g Braking index for the Vela pulsar was not determined from a standard timing analysis due 
to the large glitches experienced by this object. Instead, measurements of v were obtained 
from assumed 'points of stability' 100 days after each glitch (see Lyne et al. 1996, , for details). 



